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Abstract 



This is a companion paper of our previous work [1] where we studied the ther- 
Q^ , modynamics and phase structure of asymptotically flat black p-branes in a cavity in 



arbitrary dimensions D in a canonical ensemble. In this work we study the thermo- 
dynamics and phase structure of the same in a grand canonical ensemble. Since the 
boundary data in two cases are different (for the grand canonical ensemble bound- 
ary potential is fixed instead of the charge as in canonical ensemble) the stability 
analysis and the phase structure in the two cases are quite different. In particular, 
^ . we find that there exists an analog of one- variable analysis as in canonical ensemble, 

which gives the same stability condition as the rather complicated known (but gen- 
eralized from black holes to the present case) two- variable analysis. When certain 
condition for the fixed potential is satisfied, the phase structure of charged black 
p-branes is in some sense similar to that of the zero charge black p-branes in canon- 
ical ensemble up to a certain temperature. The new feature in the present case is 
that above this temperature, unlike the zero-charge case, the stable brane phase no 
longer exists and 'hot flat space' is the stable phase here. In the grand canonical 
ensemble there is an analog of Hawking-Page transition, even for the charged black 
p-brane, as opposed to the canonical ensemble. Our study applies to non-dilatonic 
as well as dilatonic black p-branes in D space-time dimensions. 
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1 Introduction 



Black holes in asymptotically AdS space have attracted a lot of attention in recent years. 
The reason is partly due to the AdS/CFT correspondence proposed by Maldacena [2] 
and its consequences [3, 4, 5]. The large black holes in AdS space are thermodynamically 
stable 4 [7] and so, the equilibrium states and the phase structure of such space-times can 
be easily studied. By AdS/CFT correspondence this in turn will provide us with the 
information about the similar states and phase structure on the CFT or gauge theory 
side [8, 9]. Black hole space-time is naturally associated with a Hawking temperature 
[10], so the gauge theory will also be at a finite temperature [11]. In particular, AdS 
black holes are well-known to undergo a Hawking-Page transition [7] and by AdS/CFT, 
this corresponds to the confinement- deconfinement phase transition [11] in SU(7V) gauge 
theory at large N. 

However, it is well-known that the phase structure just mentioned is not unique to 
the AdS black holes, but similar structure also arises in suitably stabilized flat as well as 
dS black holes [12, 13]. Higher dimensional theories like string or M-theory are known to 
admit higher dimensional black objects in the form of black p-branes [14, 15, 16] which are 
asymptotically flat and so, it would be interesting to see what kind of equilibria and phase 
structure they give rise to. With this motivation we obtained the equilibrium states and 
the phase structure of the black p-branes in the canonical ensemble in the previous work 
[1]. In this paper we study the same for the black p-branes in grand canonical ensemble. 
Since the boundary data in these two ensembles are quite different, their phase structures 
are also expected to be different and it is of interest to see how they differ with the change 
of the boundary data. 

We will consider the black p-branes (in D dimensions) which are the solutions of D- 
dimensional gravity coupled to a dilaton (a scalar field) and a (p + l)-form gauge field 
[15]. These are asymptotically flat and so are thermodynamically unstable 5 as an isolated 
black p-brane would radiate energy in the form of Hawking radiation [10]. In order to 
restore thermodynamic stability so that equilibrium thermodynamics and phase structure 
can be studied we must consider ensembles that include not only the branes but also their 
environment. As self-gravitating systems are spatially inhomogeneous, any specification 
of such ensembles requires not just thermodynamic quantities of interest but the place at 
which they take the specified values. In other words, we place the black brane in a cavity 
a la York [19] (see also [20, 21, 22, 23, 24]) and its extension in the charged case. The wall 

4 This is not the case for small black holes, for example, see [6]. 

5 A detailed analysis of this kind of instability for neutral black holcs/branes/rings is given in [17, 18]. 
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of the cavity is fixed at a radial distance p#, and we will keep the temperature and the 
gauge potential at the wall of the cavity (at ps) fixed. This will define a grand canonical 
ensemble (Note that for the canonical ensemble, the charge enclosed in the cavity is fixed 
instead.). So, for the grand canonical ensemble charge can vary. We will study the phase 
structure of black p-branes in this ensemble. Charged black holes in the grand canonical 
ensemble have been studied in [25, 26]. 

We will employ the Euclidean action formalism [27, 19, 28] to the dilatonic black p- 
brane geometry in D space-time dimensions. We first write the Euclidean action contain- 
ing the gravitational part including a Gibbons-Hawking boundary term [27], the dilatonic 
part and the form-field part. By using the equation of motion(i.e., on-shell), we can 
express the usual Einstein-Hilbert gravitational action plus the dilatonic part in terms 
of the form-field, therefore simplify the action and then evaluate it for the given black 
p-brane geometry placed in a cavity with the potential and the temperature fixed on its 
wall. To the leading order this action is related to the grand potential or Gibbs free 
energy of the system and is an essential entity for the stability analysis. It can explicitly 
be expressed in terms of parameters of the branes and the cavity. The stationary point of 
this action with respect to the relevant variables such as the horizon size will determine 
the relevant thermal states with the temperature and potential fixed at the wall of the 
cavity. The second derivative or derivatives of the action at the stationary point can be 
used to give us the information about the stability of the black brane at that point. We 
find the condition which determines at least the local stability of the black brane phase. 
In the case when there is a stable black brane phase, we find that there exists a minimum 
temperature below which there is only 'hot flat space' phase and no black brane phase. 
But above this temperature there exist two black brane phases with two different radii. 
The smaller one is unstable which corresponds to the maximum of the free energy and 
the larger one is locally stable and corresponds to the local minimum of the free energy. 
This local minimum becomes a global one only after the temperature rises above a certain 
value. Below this value and above minimum temperature, the locally stable black brane 
eventually makes a transition to the 'hot flat space' by a phase transition. But above this 
value, the large black brane is globally stable. Upto this point the picture is very similar 
to the chargeless black brane case in the canonical ensemble [1]. But now for the grand 
canonical ensemble as the temperature rises more, the globally stable black brane phase 
disappears after a certain value and at this point there is only one unstable black brane 
phase and the stable phase here would be the 'hot flat space'. Also note that unlike in 
the case of canonical ensemble 6 [8, 9, 12, 13, 1], here we do not have a situation where we 

6 For the case of canonical ensemble with fixed non-zero charge, the phase structure of the asymptoti- 
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have two locally stable black brane phases and so there is no such phase transition similar 
to the van der Waals-Maxwell liquid-gas phase transition or even a transition from one 
black brane phase to the other. We would like to remark that in this grand canonical 
ensemble, unlike in the case of AdS black holes [8, 9], we do not have a critical value 
for the potential in the sense defined there since the existence of this critical value in 
AdS case is largely due to the presence of the cosmological constant ( as an additional 
parameter). One can check easily that if we set the cosmological constant vanish (i.e., 
setting the parameter / — > oo there), the extremal black holes becomes supersymmetric 
BPS ones, which is in analog of the present discussion, and the potential will not have the 
critical value anymore but has only an upper bound. The existence of the upper bound 
of potential for the present case is explained below and in section 3. 

We now discuss the case of extremal black p-branes. Mathematically, the extremal 
brane configuration can be obtained simply by taking the 'mass-equal-charge' limit of the 
corresponding non-extremal one. However, this does not mean that the extremal branes 
and non-extremal branes can actually be related to each other by a physical process which 
can also be understood from the semi-classical discussion given by Gibbons and Kallosh 
[29], Hawking et al [30, 31] and Teitelboim [32], for the case of black holes including 
the dilatonic ones (the latter resembles the present discussion). So, for example, the 
entropy changes discontinuously from the non-extremal case to the extremal one. The 
non-extremal case has a non- vanishing entropy, given always by one quarter of the event 
horizon area while the extremal one always has vanishing entropy 7 . This very fact, just 

cally flat black p-branes [1], when d > 2, with d = D — 3 — p, resembles the one in [8, 9] for charged AdS 
black holes as discussed in detail in [1]. For examples, in both cases, there exist a critical charge and when 
the charge carried by the branes is greater than this charge, there always exists a thcrmodynamically 
stable state for each given temperature. When the charge is less than the critical value, there exists a line 
of first order phase transition of van der Waals-Maxwell liquid-gas type, ending at a second order phase 
transition point when the charge reaches the critical value. The critical exponent for the specific heat is 
universal and has the value —2/3. However, for d = 1 the phase structure is similar to the chargeless case 
except that there does not appear to exist an analog of Hawking-Page transition due to the presence of 
charge and the horizon size has a lower bound set by the charge (rather than zero as in the chargeless 
case). The d = 2 case is a borderline between d = 1 case and d > 2 cases. Similar to the cases of d > 2, it 
has a critical charge and three subcases depending on whether its charge is greater than, equal to or less 
than the critical charge. However, there exists neither a critical behavior nor a Van der Waals-Maxwell 
type liquid-gas phase transition. The stable phase structure looks more like the d = 1 case. For other 
details we refer the reader to [1]. 

7 The discrepancy in the entropy between the semi-classical calculation and the microstate counting in 
string theory in the AdS/CFT context has been addressed recently in [33]. It was pointed out there that 
the non-vanishing of the entropy is due to a different reason and the semi-classical result of vanishing 
entropy for the extremal black hole still holds. 
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like the case of black holes, implies that the non-extremal branes and extremal branes 
are qualitatively different objects and a non-extremal black brane cannot turn into an 
extremal brane. Actually, the nearer the mass of the non-extremal brane gets to the charge 
the lower the temperature of the brane and so the lower the rate of radiation of mass. Thus 
the mass will never be exactly equal to the charge. This, in turn, implies, as it should, 
that there does not exist any real physical process which can bring down the temperature 
of a thermal system to absolute zero since otherwise it would be against the third law 
of thermodynamics that the zero temperature can never be reached in reality. Moreover, 
also like the case of black holes, an extremal brane cannot become a non-extremal brane 
even when matter or radiation are thrown into the extremal one, which seems at first 
sight contrary to the common sense [30]. This is due to the fact that one can identify 
extremal branes with any period and therefore they can be in equilibrium with thermal 
radiation at any temperature. Thus they must be able to radiate at any rate, unlike non- 
extremal branes, which can radiate only at the rate corresponding to their temperature. 
So extremal objects always radiate in such a way as to keep themselves extreme when 
matter or radiation is sent into them. This discussion applies to the asymptotically flat 
branes in canonical ensemble for which the charge inside the cavity is fixed. 

In this paper, we are considering the asymptotically flat branes in grand canonical en- 
semble for which the potential at the wall of the cavity is fixed instead. Before we discuss 
the thermodynamical stability, we first need to address the question: given the bound- 
ary data in this ensemble mentioned earlier, what are the possible stable configurations 
allowed in the cavity which can be used to discuss this stability? Unlike in the canonical 
ensemble, the charge inside the cavity in the grand canonical ensemble is not fixed and can 
fluctuate or exchange with the environment. So one expects first that the 'hot flat space' 
can be a stable phase under certain conditions. The second candidate is the non-extremal 
branes, charged or uncharged, and the third possible one is the extremal branes. However, 
the probability for the existence of non-extremal branes is much higher than that of the 
extremal branes for the following reasons. In the case of black holes, the pair creation rate 
of the non-extremal black holes is enhanced by a factor of e s over that of the extremal 
branes where S is the entropy of the non-extremal black holes [30, 29]. We expect this 
to apply to the case of branes as well. Further, unlike a non-extremal brane, an extremal 
brane with a large charge can fission (split) into extremal branes with smaller charges 
without violating the second law of thermodynamics since this process is suppressed by a 
factor of e AS with AS the change of total entropy of the system of extremal branes and 
for extremal branes AS = 0. Note that an extremal brane with a macroscopic charge 
in the present context (for the gravity description to remain valid) is actually a BPS 
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brane, preserving one half of the underlying spacetime supersymmetries and there is no 
binding energy among BPS branes with smaller charges, even the individual BPS branes 
with unit charge (which constitutes the extremal branes with macroscopic charge). So a 
large extremal branes can easily split into smaller extremal ones. This provides another 
means to understand the above fission process. In other words, the extremal branes are 
unstable in the first place in this ensemble and should be excluded from consideration in 
our following thermodynamical stability analysis 8 . 

This paper is organized as follows. In section 2, we discuss the dilatonic black p-brane 
solution and evaluate the action in Euclidean signature. Then we discuss the stability 
of various equilibrium states from this action in section 3. The phase structure of the 
black p-brane is discussed in section 4. Then we conclude in section 5. A more general 
two- variable stability analysis (as opposed to one-variable stability analysis performed in 
section 3) is presented in the Appendix. 



2 Black £>-brane solution and the action 



The black p-brane solution was originally constructed [14] as a solution to the ten di- 
mensional supergravity containing a metric, a dilaton and a (p + l)-form gauge field and 
was generalized to arbitrary dimensions in [15]. These solutions are given in Lorentzian 
signature, but for the purpose of studying thermodynamics, we write the black p-brane 
solution in Euclidean signature as (see for example [34]), 



d-l 



a + a_ d - 2 dt 2 + a°~ 2 Y^( dxi f + a;'a m dp 2 + P 2 Ai d dn 2 s+1 



i=l 



A 



a0 o /2 



J/2 



r_r i 



d/2 



dt A dx 1 A ... A dx p , 



F [p+2] 
X<t>-<M 



dA [p+1] 



-ie a * o/2 d ( r - T - +y * d/ ~ dp Adt A dx 1 A ... A 



dx p , 



Here we have defined 



P 



(2) 



8 So the natural background in the grand canonical ensemble in the present discussion is the 'hot 
flat space' while from the previous discussion the natural background in the canonical ensemble is the 
extremal branes. 
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where, r± are the two parameters characterizing the solution and are related to the mass 
and the charge of the black brane. In the metric (1) the Euclidean time is periodic and 
so, the metric has an isometry S 1 x SO(d— 1) x SO(d + 2) indicating that it represents a 
(d — 1) = p-brane in Euclidean signature. The total space-time dimension is D = d+d+2, 
where the space transverse to the p-brane has the dimensionality d + 2. <fi is the dilaton 
and (f)Q is its asymptotic value and related to the string coupling as g s = e*°. a is the 
dilaton coupling and is given for the supergravity theory with maximal supersymmetry 
by, 

a 2 = 4 . 3 

D-2 V ; 

It is clear from the Lorentzian form of the above metric that when r_ = 0, and a = 0, it 

reduces to the D- dimensional Schwarzschild solution which has an event horizon at p = r + , 

whereas, at p = r_, there is a curvature singularity. So, the metric in (1) represents a 

black p-brane only for r + > r_, with r + = r_, being its extremal limit [14]. A p-brane 

naturally couples to the (p+l)-form gauge field whose form and its field strength are given 

in (1). Note that we have defined the gauge potential with a constant shift, following [26], 

in such a way that it vanishes on the horizon so that it is well-defined on the local inertial 

frame. The black p-brane will be placed in a cavity with its wall at p = ps- It is clear 

from the metric in (1) that the physical radius of the cavity is 

p B = A?p B (4) 

while ps is merely the coordinate radius. So, it is pB which we should fix in the following 
discussion and not ps- Note that when the black brane is non-dilatonic a = and in 
that case ps = Pb- Also we fix the dilaton 9 at ps, which indicates that the asymptotic 
value of the dilaton is not fixed for the present consideration and this is crucial for our 
discussion in expressing relevant quantities in 'barred' parameters. By this argument we 
also have 

f± = A«r ± (5) 

and r± are the proper parameters which we should use in the present context. In terms 
of the 'barred' parameters A± remain the same as before, 

A ± = 1 - % = 1 - !4. (6) 

PB PB 



9 This enables us to obtain its correct equation of motion from the corresponding action in the presence 
of a boundary. 
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Since the Euclidean time coordinate in (1) is periodic so, for the metric to be well-defined 
without a conical singularity at p = r + , the Euclidean time must have a periodicity, 

which is the inverse temperature at p = oo. The local (3(pb) is given 10 as 



i 



(3 = f3(p B ) = A|A_ 2(d+d, /3* (8) 

which is the inverse of local temperature at ps when in thermal equilibrium with the 
environment (the wall of the cavity) at the inverse temperature (3. Note that in terms of 
the 'barred' parameters the inverse of local temperature P(pb) can be expressed from (8) 
and (4) as, 



ii . .ii 



where in the second equality A± are also expressed in terms of 'barred' parameters. Also 
the charge is defined as, 



Qd = 4- /e- aW *F b+2] = %e-^/ 2 J(r + r_)^ 
V2k J V2ft 

' d e- a ^\f + f.f/ 2 . (10) 



V2k 

In (10), k is a constant with l/(2n 2 ) appearing in front of the Hilbert-Einstein action in 
canonical frame but containing no string coupling g s , *F[ p+ 2] denotes the Hodge dual of 
the (p + 2)-form field given in (1). Also, Q n denotes the volume of a unit n-sphere. In 
the last line of (10) we have expressed the asymptotic value of the dilaton in terms of the 
fixed dilaton via the boundary condition <P{Pb) = at the wall of the cavity from the 
relation (1) and then expressed r± by f± from (5). 

In the grand canonical ensemble, the fixed quantities are the physical radius of the 
wall of the cavity pb, the temperature T, the brane volume V p (The local volume at ps, 

d(d-l) 

Vp(Pb) = A^ (D_2) V*, with V* = f d p x, is set equal to the fixed value V p ) and the potential 
$, all at the wall of the cavity. They are independent parameters and the corresponding 
local quantities for the branes have to be equal to these pre-selected fixed values on-shell 



10 Thc corresponding local temperature is related to the temperature at spatial infinity by the so-called 
Tolman relation [35]. 
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at the wall of the cavity, respectively. The potential in the local inertial frame at the wall 
of the cavity can be obtained from ^4[ p +i] and the metric in (1) as, 

A [p+1] = -ie a ^ 2 (A_A + )~5 2 U _ 1+ J dth dx 1 ... A dx p 

= -iV2Mi A dx 1 . . . dx p (11) 
where (i, x 1 , . . . , x p ) are the coordinates in the local inertial frame and is related to the 



original coordinates as i = Ajj:A_ 2(D ~ 2) t and x % = A^ D ~ 2) x l for % = 1, 2, ... ,p as can be 
seen from the metric in (1). So, $ is the potential conjugate to the charge and is set equal 
to the fixed value <3> at the wall. We will now evaluate the action with these boundary 
data. 

The Euclidean action for the dilatonic black branes in the canonical ensemble has 
already been evaluated in the Appendix of our previous work [1]. We will use that result 
to evaluate the action for the grand canonical ensemble. The relevant action for the 
gravity coupled to the dilaton and a (p + l)-form gauge field in the canonical ensemble is, 

I% = I%(9) + !%(<!>) + (12) 

where, I gig) is the gravitational part of the action which has a Hilbert-Einstein term 
and a Gibbons-Hawking boundary term [27], Ie{</>) * s ^ ne dilatonic part and Ie(F) is 
the form- field part. The first two terms remain exactly the same in the grand canonical 
ensemble, however, the form-field part which in the canonical ensemble has the form, 

j£(F) = -!___L_ / d D xJge- a ^F 2 d+1 

I ^^VT^e-^F" 2 -^,,,.,, (13) 

*K « ! JdM 

will differ in the grand canonical ensemble since the potential at the wall of the cavity 
is fixed and so, the last term in (13) will be absent in the grand canonical ensemble. 
Therefore, the action for the grand canonical ensemble can be obtained from that of the 
canonical ensemble by the relation, 

I e° = ^ + ^ [ d D ^x^n,e- a ^F^-^A^ 2 ..., d . (14) 
zk ai j dM 

Here dM denotes the boundary of the space-time M and n M is a space-like vector normal 
to the boundary. 7 is the determinant of the boundary metric. Note that the spacetime 
M here is the bulk region defined by ps > p > f+ while the boundary dM is defined as 
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the surface at p = pB plus the one at p = f + . The surface integration in (14) at p = f+ 
vanishes due to the vanishing form-potential there. Evaluating the last term in (14) for 
the black p-brane configuration given in (1) we get, 

1% C = I°e- ^n s+1 drlA + = Jg - f3V p Q d ^>, (15) 
where we have used the charge expression and the potential from (10) and (11) as, 

Q d = %^e-^ 2 (f + f_)i, 



y/2. 



K 



i 



V2k V r +y V A - 

Now using the explicit form 11 of the Euclidean action for the canonical ensemble given in 
[1] we write (15) as 12 

I E (P,$,p B ,V p ;Q d ,f + ) = f3Q 

= pE(p B , V p , Q d , f+) - S(p B , V p ; Q d , f+) - (3V p Q d $> 

" ~^2^ Pb 



/A \ 1/2 

(d + 2) f^±J +J(A + A_) 1 / 2 -2(J+l) 




2k 2 

where Q = E — TS — Q§ is the Gibbs free energy. In the above, we have the total charge 
enclosed in the cavity as Q = V p Q d , the entropy 

2k 2 + ~ 




11 In obtaining this, we have made use of the on-shell equation of motion 

R _ I( 5 0) 2 = ^— - - e -^F\ (17) 

2 V v ' 2(D-2)(d+l)! V ; 

to express Ie{!j) + Ie(<P) m terms of the form-field F to simplify the computation. 

12 In obtaining the explicit action, we have used the background subtraction approach. The natural 
background here is taken as the 'hot flat space' with a constant dilaton (j>, a fixed potential $ and other 
boundary data. The equations of motion are satisfied trivially inside the cavity for the 'hot flat space' 
with the given boundary data. 
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where V* = J d p x as defined before and is independent of the location of the cavity as it 
should be, and the energy for the cavity as 



E 



Pb 



2k 2 



(d + 2) 



A_ 



1/2 



d(A+A. 



,1/2 



2(d+l) 



(20) 



which gives the ADM mass when ps — >■ oo. Since this is the Euclidean action for the 
black p-brane in the grand canonical ensemble we will be working with, for brevity, we 
have removed the superscript t GC from Ie in writing (18). Note that we have expressed 
everything on the r.h.s. in terms of the 'barred' parameters showing that the formalism 
works for both the non-dilatonic as well as dilatonic branes. In (18) /3, Pb, $ are the 
inverse of temperature, the physical radius and the potential of the cavity, therefore are 
all fixed, and Qd, f+ are variables. Note that f_ is not independent and can be expressed 
in terms of Qd and f+ from (16) as, 



V^Qd a} 



(21) 



Substituting (21) in (18) we can express the action in terms of two variables Qd and f+. 
Then varying the action with respect to Qd and putting that to zero, i.e., at the stationary 
point we obtain (after some simplification), 



diE 

dQ 







alL d+l' + V 



A + V 



d I 2 . i 
d 



\ 



4vrf + AiA d 

1 

ii x 



2 d 



(22) 



Similarly, varying the action with respect to the other variable r + and setting that to 
zero, i.e., at the stationary point we obtain, 



dip Anr i i — 4 / r 

|^ = => fi = p{p B ) = ^±A| A_^ [ 1 - ^ 



d \ 3 2 



(23) 



This is the correct form of (3 we had given earlier in (9). Using this (23) in (22) and 
substituting the form of Qd given in the first equation of (16) we recover the correct form 
of the potential given in the second equation of (16). This is a verification that the action 
Ie given in (18) is indeed correct. We will use this form of the action in the next section to 
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study the stability of the equilibrium states of the black p-branes in the grand canonical 
ensemble. 

Before we close this section, we would like to discuss briefly the validity of the use 
of the effective action in describing the phase structure of black p-branes throughout the 
parameter space in the present ensemble. Since we are considering a situation where the 
black p-branes are placed in a cavity, this implies that the size of the cavity pB must be 
larger than the horizon size f+ of the branes. In obtaining the explicit Euclidean action 
(18), we have performed the integration of the bulk action in the range of ps > p > r+- 
To have a valid description in terms of the effective action, we need to keep the curvature 
of black brane spacetime uniformly weak throughout this range and if the fundamental 
theory is string theory, we need in addition to keep the effective string coupling uniformly 
weak. One can check that the requirement (s) can be satisfied with the expected condition 

f+ = r + A a y Ad » /, (24) 

where I is the fundamental length scale of the underlying theory, for example, it is the 
Planck scale l p in eleven dimensions or the string scale l s in ten dimensions. 

There appears no further subtlety in the grand canonical ensemble considered in this 
paper since the extremal branes are excluded from consideration as discussed in the In- 
troduction. The relevant configurations here are the 'hot flat space' and the non-extremal 
branes. Note also that the location of the cavity wall along the radial direction satisfies 
pB > p > f+ > f_. However, there appears to have a potential issue regarding the 
fixed dilaton value at the wall of the cavity in the canonical ensemble, considered in our 
previous work [1], for certain extremal branes such as D3 branes for which the so-called 
attractor mechanism is at work 13 . Once this mechanism is at work, the dilaton value at 
the horizon is completely determined by the relevant charge(s), independent of its asymp- 
totic value. The extremal p-brane configuration from the supergravities with maximal 
SUSY in spacetime dimension D > 4 can be obtained from the non-extremal one given 
in (1) by taking the extremal limit, i.e., r + = r_. These branes preserve one half of 
the space-time supersymmetries and the dilaton at a given p is related to its asymptotic 
value via the relation given in (1). The (arbitrary) fixed value at the wall of the cavity is 
therefore related to the (arbitrary) asymptotic value of the dilaton with the fixed cavity 
size. So in the present context, one can re-phrase the attractor mechanism as: the value 
of the dilaton at the horizon is independent of the fixed value of the dilaton at the wall 
of the cavity. So if the attractor mechanism is at work, the wall of the cavity cannot 
be too close to the horizon since otherwise the approach will be inconsistent with this 

13 We thank the anonymous referee for raising this point to us. 
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mechanism or it will give rise to a jump of the value of the dilaton in the bulk which 
is impossible. Fortunately, the situation is not so restrictive as indicated above. Once 
the attractor mechanism is at work, any given point outside the horizon has actually an 
infinite physical distance away from the horizon (see, for example, [36, 37, 38]) and so the 
condition ps > r + = r_ with fixed pb, the radial location of the wall of the cavity, will 
be enough. 



3 Stability analysis of the black £>-branes 

For the purpose of the stability analysis we will define some new parameters following 
refs.[26, 13, 12] and rewrite the action, the potential and the inverse temperature in terms 

of those parameters. Let us first define the reduced charge as, 

i 



Then from (21) we have 



r+ 

Next we define the following parameters, 



(26) 



I = (£tV, i = ? (27) 

\Pb J 4npB \Pb J 

where the dimensionless parameter b is fixed, but the other two parameters x and q vary. 
Note that the parameter b is related to the inverse of temperature of the environment, q 
is related to the charge and x is related to the horizon size. In terms of these parameters 
we have 

fid, 

A+ = 1 - -t = 1 - x 

P% 

A 1 _^1 = 1 _ ( W~ _ 1 <f 



Pb KPb x 



r d 2 

Using (28) the two equations at the equilibrium (at the stationary point) giving the 
constant temperature and the potential at the wall of the cavity given in (23) and (22) 
can be written as 

b = b(x,q), (p = (p(x,q) (29) 



13 



where we have defined </? = \f2ne a< ^/ 2 <3> and 

b(x, q) = \- X ^ , X - T , <p(x, q) = - I 7 ^ I 




r/ 



The two equations in (29) can also be derived from the following reduced action in a 
similar fashion as (23) and (22) but now with respect to variables x and q, respectively, 



I E (x,q) = 



2k'I 



47ip d B +1 V p Q 




> (31) 



where Ie is as given in (18). These two equations, for given b and <p, determine x and 
q completely. However, in the presence of two variables, the analysis of stability at the 
extremal points determined by these two equations is a bit more complicated than that 
given in [1], in the canonical ensemble. There the relation at the charge equilibrium 
for fixed charge has been employed to reduce the two variables to only one 14 . While 
the stability analysis with two variables, which we will perform in the Appendix, is more 
complete, there actually exists an analogous simpler analysis as in the canonical case. But 
for this we need to employ the potential relation at the equilibrium for fixed potential, 
as given in the second equation of (29), to reduce the two variables to one variable. In 
what follows, we will first perform the simpler one- variable analysis and find the stability 
condition. We will perform the more complete yet more complicated two-variable analysis 
in the Appendix and show that these two approaches give the same stability condition. 

Note that since x/q = (f + /Q^) d = (f + /f„) d / 2 > 1, so, x > q and since x = (f + /pB) d < 
1, so, a: has the range q < x < 1. The condition (24) may put a more strict lower bound for 
x if q is too small. The end point x — 1 corresponds to taking pB — f + , i.e., the cavity is 
placed on the horizon where the ensemble temperature is infinity. In practice, we can keep 



14 While the one- variable stability analysis given in [1] appears natural in the canonical ensemble, there 
exists also a more general two- variable analysis which is more complete yet more complicated and gives 
the same stability condition. For the present case, the two- variable analysis seems more natural but as 
we will show in the text that when the second equation in (29) for the fixed potential is used to reduce 
the two variables to only one, the same stability condition can be obtained. 
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an ensemble at any given high temperature and the infinitely high temperature should 
be understood as a limiting process. So we should take x < 1. The other end point 
x = q corresponds to the extremal case which should be excluded from the consideration 
given what has been discussed in the Introduction. For this reason, we take x > q. So, 
following [26], we limit ourselves to the physical region ofO<g<a;<lin what follows, 
keeping in mind the condition (24) 15 . Also note from tp expression in (30) that since 
1 — q 2 /x = 1 — (q 2 /x 2 )x > 1 — x, so we have < <p < 1, with ip = corresponding to 
the chargeless (q = 0) case. As discussed in Introduction, we should exclude the ip = 1 
(x = q) since it corresponds to the extremal case, ip > 1 cannot be in a black brane phase 
and it actually corresponds to the 'hot flat space' phase. 

Let us now consider the one-variable case first, i.e., when the second equation in (29) 
< <p = <p(x, q) < 1 is satisfied. We can now solve this equation to get, 

q - = ^ (32) 
x 2 1 - (1 -<p 2 )x 1 ' 

Now substituting (32) in b(x, q) equation in (30) we have for fixed potential <p, 

b{x) = b(x, q) = : — , (33) 

d(l -<^ 2 )2"3 

and in the reduced action in (31), we have 
I E {x) = I E (x, q) = -2b (d+ l) [(1 - (1 - f) x) 1/2 - l] - x 1+1 ^ (1 - ^ 2 ) 1/2+1/J . (34) 



One can easily check that 

dI E {x) (l+rf)(l-<^ 2 ) 



[b-b(x)], (35) 



dx [1 - (l -<f?)x] 1/2 
where b(x) is defined in (33). This vanishes at the stationary point of the action, i.e., 

b = b{x), (36) 

which is the first equation in (29) in the present consideration. The local stability at the 
stationary point is determined by requiring it to be a local minimum of the action, i.e., 



d 2 I F (x) 



dx 2 



d)(l-^ 2 ) db{x) >Q 



[1 - (1 -<p 2 )x} 1/2 dx 
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As will become clear, the local stability condition given in (41) does satisfy the condition (24). 
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In other words, we need to have (since the other factors are all positive) 

2- (2 + d)x(l -tp 2 ) 



db(x) x ^ 



< 0, (38) 



dx 2d 2 {l - £ 2 )V2-i/d [1 - (1 - v9 2 )x] 1/2 
which gives the local stability condition as 

2 - (2 + d)x(l - <^ 2 ) <0. (39) 

In the above, x is a solution of the present equation of state as given in (36). 

As we will show in detail in the Appendix that the above local stability condition (39) 
continues to hold when we consider the more general yet more complicated two-variable 
situation and this justifies the simplified one- variable analysis of stability performed above. 

We therefore conclude that for the local stability of the system we must have (39) to 
be satisfied. To understand the meaning of the stability condition (39) we rewrite it as 

x (1 - f) > - J— (40) 
K J d + 2 

which further implies, 



< x < 1, and < ip < \\ ~ - . (41) 
d+2 V d+2 

In other words, for given (p satisfying the above constraint, only x which is in the range 
specified in (41), gives at least a locally stable system. Otherwise the system is not 
stable. Note that the lower bound 2j{d + 2) for x is of the order of unity for allowed 
d < 7 and the condition (24) in the present context is 1/pB 1. So we can easily have 
1 > x > 2 /(d + 2) > 1/pB to satisfy the condition (24). 

This concludes our discussion on the stability of asymptotically flat black p-branes 
in the grand canonical ensemble. Using this information we will construct the phase 
structure of the equilibrium states of the black p-branes in this ensemble and compare 
with that in the canonical ensemble. 



4 Phase structure of black £>-branes 

We have seen in the previous section that the stability condition for the black p-branes 
is given by (39). When we use the potential condition at the equilibrium, (p = (p(x, q), to 
eliminate q and obtain (33), i.e, 

d(l -02)3-3 
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we have the range of x as < x < 1. This is different from the fixed charge case of the 
canonical ensemble where the range is q < x < 1, since here we are considering the grand 
canonical ensemble where the potential (p is fixed and not the charge q. Even for the zero 
charge case of the canonical ensemble, where the range of x in the two cases are the same, 
the phase structure in these two cases, as we will see, will be quite different due to the 
different boundary data in the two ensembles. Note that as x — > 0, b(x —>•())—>■ 6(0) = 0, 
exactly as in the canonical case, but as x — > 1, b(x — > 1) — >■ 6(1), where it has the value, 



6(1) 



i i 



> 



(43) 



d(l-(p 

and this is different from the fixed charge canonical case including the zero charge case 
The local temperature at this point has the value 

d(l 



T(l) 



i_ i 



(44) 



47T p B <f 

Also note from (42) that in the range < x < 1, 6(x) > and so, if we plot b(x) vs. x 
curve it will start from zero when x — > and then it rises and ends at 6(1) > given 
above when x — > 1. In between there can be extrema for b(x) and in fact, it is not difficult 
to check from 



db(x) 



b(x) 2 - (d + 2)x(l - if) 



(45) 

dx 2dx[l -x(l -ip 2 )} 

that there can be only one extremum and that is a maximum if it exists at all. This max- 
imum is not always guaranteed to exist unlike the chargeless case of canonical ensemble. 
However, when it exists indeed, the characteristic behavior of b(x) vs. x is given in Figure 
1. The plot of T(x) vs. x curve will have opposite behavior. It will start from infinity 
when x — > then the curve will drop and end at T(l) as x — > 1. In between there can be 
an extremum and in this case it is a minimum. Let us now determine the condition for 
the existence of the maximum in b(x) or minimum of T(x). Equating (45) to zero gives 
us the maximum of x as, 

(46) 



X i ■ 



(?) 



(2 + d)(l 

Now substituting this value in (42) we determine the maximum value of b(x) and from 
there the minimum value of T(x) as, 



T 




d(d + 2) (l-f) 
d(d + 2) (1 



(47) 
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Figure 1: The typical behavior of the (reduced) inverse temperature b(x) vs. the (reduced) 
horizon radius x at a fixed potential (p. 

The condition for the existence of maximum is x max < 1, which gives <p < (d/(d + 2)) 1 / 2 , 
i.e., the same constraint we obtained earlier in (41), which is expected since the non- 
existence of maximum implies no local stability. From this we have the following constraint 
on 6 max or T min , 

1 _ i 

1 ( 2 \ 3 m V2d /2 + dV 

6 max < —= = T min > — — -— 48 

V2d\2 + dJ 47rp B ^ 2 J 

where we have used the constraint on (p given in (41). So, when 6 max (T min ) given in 
eq.(47) satisfy the constraint (48) we will have the maximum (minimum). From Figure 1, 
we can see that when x max < x < 1, db/dx < and we have a locally stable system and 
when < x < x max , db/dx > and we have an unstable system consistent with what we 
discussed earlier. 

For any given (p < 1 and b, satisfying 6(1) <b< b max with 6(1) as given in (43), one 
expects two black brane solutions with radii x\ and X2, where x\ < X2 from b = b(x) as 
shown in Figure 1 and only the large one (£2) will give a locally stable phase. The small 
one (xi) will be unstable. Since here for each given (p and b, we have only one locally 
stable phase, we don't have a thermodynamical phase transition like in the canonical 
ensemble 16 , even though x\ = 22 = x max seems to appear as a critical point. When 
<P > (d/(d + 2)) 1 / 2 and/or b > 6 max (T < T min ) no black brane phase is possible in the 

16 In the canonical ensemble there is a phase transition between two stable black brane phases, one 
of which is locally stable and the other is globally stable. Which one is locally stable and which one is 
globally stable depend on the temperature of the environment in contact with the system [1]. 
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grand canonical ensemble and the only possible thermally stable phase here is the 'hot 
flat space'. Also for = 6(0) < 6 < 6(1), assuming <p < (d/(d + 2)) 1//2 from now on, no 
stable black brane phase is possible and the 'hot flat space' is again the thermally stable 
phase. This clearly shows that the thermodynamics and the phase structure are quite 
different for the grand canonical ensemble than those in the canonical ensemble [1] and 
this is expected since the boundary data are different in the two cases. 

To the leading order, the Euclidean action is directly related to the grand potential or 
Gibbs free energy as f3Q = Ie- So, for a given temperature, smaller the value of Ie, the 
more stable the system is. We will consider the reduced action Ie- For a given 6 and (p 
satisfying their respective constraints such that the system is locally stable, if the reduced 
action is positive at the stationary point, it is metastable since the reduced action has 
smaller value, i.e., zero, for the 'hot flat space' with the same boundary data and will 
make a transition to this phase. In other words, if the system is initially at the locally 
stable black brane phase, after some time, a phase transition would occur so the black 
brane phase will become the 'hot flat space' phase with a different topology via a phase 
transition. So, for making sure that we have a stable black brane phase, we need to look 
for condition such that the minimum of the reduced action is actually a global minimum. 
This can be simply realized by conditions such that the stationary action at the minimum 
is negative. 

The reduced action at the stationary point can be expressed from the action (34) as, 

lE= b -(d + 2)(l-y)(y-J-) (49) 

y V a + 2 J 

where we have used b(x) = 6 with x the stationary point and y = Wl — x(l — <p 2 ). Given 
the fact that y < 1, so a negative reduced action requires, 

y < J- (50) 
d + 2 

which in turn gives, 

x(l-^)> A -i±^, (51) 
K J (d + 2) 2 K J 

where we have used the expression of y given above. Note that this condition is consistent 

with the local stability condition given in (40) since A(d + l)/(d + 2) 2 > 2/(d + 2). Also 

note that 4(d + 1)/ (d + 2) 2 < 1 for all d and therefore, there is no contradiction. Actually, 

the constraint (51) gives the following new constraints 

4(J+1) d , . 

-X '- < x < 1, 0<w< ^ . 52 

{d + 2) 2 d + 2 
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So (51) or (52) is the condition for global stability while (40) or (41) is merely the local 
stability condition. If we denote 



* B ~ (d+2ni- ^) {> w) (53) 

then for < <p < d/(d + 2) and 6(1) < 6 < b g (T g < f < T(l)), where, 



(rf + 2) 1+ 3 v /f- <p< 



Fx 

and 



4(d+l) 

b g = - v a / < 6 max (54) 



(d+2) i+ i^r^2 ^ 

— ~ -'mill) \<- ) <- ) J 



Airp B U(d 



the reduced action has a global minimum at x g < x^ < 1- In other words, this ensemble 
can give a sensible description of black brane thermodynamics only when the above con- 
ditions are satisfied. We would like to point out that when the maximum of b(x) exists 
(or minimum of T(x) exists), the grand canonical ensemble is in some sense similar to the 
chargeless case of canonical ensemble except around x — 1. 

While the above gives a clear picture of phase structure for black p-branes in the 
grand canonical ensemble, we can also consider the dependence of the free energy on 
the temperature at fixed potential, as in the usual approach, to understand the same 
structure. For this, we note the known results from the above: 1) When (p > 1, the only 
stable phase is the 'hot flat space' and we have the vanishing free energy (Ie = 0). 2) 
When (d/(d + 2)) 1 / 2 < (p < 1, we have the (reduced) on-shell free energy for the branes 

fi(, )s ^ = ( J +2) iz»f v _ ^L) >0 , (56) 
b y \ d+2 



where y = ^1 — x(l — tp 2 ) and (d/(d + 2)) 1 / 2 < y < 1. When this happens, b(x) does 
not have a maximum and there exists only one solution x from the equation b = b(x) for 
= 6(0) <b< 6(1). The reduced temperature defined as t = 1/6 = AnpsT is now in the 
range of t(l) < t < oo with 

from (43), and the plot of the (reduced) free energy vs the (reduced) temperature gives 
only the unstable branch as shown in Figure 2(a). For this case, we do not have stable 
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brane phase and the only stable phase is the 'hot flat space'. Here we have Cl(t — » oo) — > 
and 



Cl(t t(l)) -> (d + 2)^— - | 



(5* 



which is positive for the allowed range (d/(d + 2)) 1 / 2 < (p < 1. Let us examine the slope 
of S7(t) with respect to t in general. We have 



dCl(t) 
dt 



1-^ 

2y3 



^2 



x(l-^ 2 ) 



d + 2 



dx(t) 
dt 



d x{t)(l-$ 2 ) 



<0, 



(59) 



d + 2 *2/(*) 
where in the second line we have used 

&(*)* = 5* = 1, ^P^l 

ax at 

and (45). In other words, the slope is always negative. 3) When < (p < (d/(d + 2)) 1 / 2 , 

b(x) has a maximum and occurs at x max given in (46). This gives a minimum (reduced) 
temperature as 



(60) 



2 + d 



d(d + 2) (1 -^ 2 ) 



(61) 



There exist two solutions x\(t) and X2{t) (x\ < X2) from b(x) = 1/t for each given t 
satisfying t(l) > t > t mm as discussed earlier and shown in Figure 1. The small xi(t) 
corresponds to the unstable brane and the large X2(t) corresponds to the locally stable 
brane. This implies that there always exist two branches in the plot of the (reduced) free 
energy vs the (reduced) temperature. The unstable branch is obtained by substituting 
the solution Xi(t) in the reduced free energy (56) with t min < t < 00 while the locally 
stable branch is obtained similarly but now by substituting the solution x%(f) instead with 
tmin < t < t(l), where t(l) is again given by (57) but (p has its present range as given 
above. The unstable branch always has positive free energy which is larger than the free 
energy in the locally stable branch except at t — t m m where the two become equal and it 
is 



^(t = tmin) = (d + 2) 1 




> 0, 



(62) 



which is explicitly independent of ip in the present case. The (reduced) free energy at the 
other end in this branch is again Q(t — > 00) — > 0, also independent of (p. In this branch, we 
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- t 



(a) (b) (c) 



Figure 2: The typical behavior of the (reduced) free energy Q(t) vs. the (reduced) tem- 
perature t at a fixed potential (p. Plot (a) considers the case \J d/(d + 2) < <p < 1 which 

has only one unstable branch. Plot (b) considers the case d/(d + 2) < (p < \J d/(d + 2) 
which has two branches, one unstable and the other only locally stable branch. Plot (c) 
considers the case 0<(p<d/(d + 2) which also has two branches, one unstable and the 
other (locally) stable branch. The unstable branch has always positive free energy for 
tmir, < t < oo while the stable branch has negative free energy in the range t g < t < t(l), 
where it becomes globally stable. 
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always have xi(t)(l — <p 2 ) < 2/ (d + 2) which can be obtained from xi{t) < x max with x max 
given in (46). This branch gives no stable brane phase. For the locally stable branch, 
when d/(d + 2) < (p < (d/(d + 2)) 1/2 , Q(t) > for t min < t < t(l) with the positive 
Q(t = t m in) given by (62) and the non-negative Q(t — > t(l)) still given by (58) but now 
for the present range of tp given above. The plot of the free energy vs. the temperature 
is given in Figure 2(b) for this subcase. This locally stable branch gives only meta-stable 
brane phase and will eventually decay to the 'hot flat space' after their formation. If 
we want the existence of globally stable branes, the free energy needs to be negative at 
certain temperature. This can be realized when < (p < d/(d + 2). With this, the free 
energy vanishes at t = t g and becomes negative when t g < t < t(l) with 

(j+2) i+ iy/r^ 

tg = l/bg = J , (63) 

(4(J+1)) 3 

from (54), therefore this branch becomes globally stable in this range. One can check 
that the (reduced) free energy Q(t — > t(l)) < as expected in the stable branch using its 
expression given again by (58). Note that the value of this limit increases as (p increases. 
One can also check that t min <t g < t(l) in the present subcase. The plot of the (reduced) 
free energy vs the (reduced) temperature for fixed <p is now given in Figure 2(c). For each 
given t satisfying t min < t < t(l) in the present two subcases, the corresponding plot just 
reflects what has been described. In the range < (p < d/(d + 2), t m i n ,t g and t(l) all 
decrease when (p increases which can be seen from their respective expressions. Note that 
&(t = tmin) i n both branches and Cl(t — > oo) in the unstable branch are both independent 
of <p, and Q(t — > t(l)) increases as <p increases. As shown above, the slope of Q(t) is 
negative in either branch. These features give the general character of the free energy vs 
the temperature plot in both branches for different fixed potentials. One may take Figure 
2(b) as the consequence of increasing <p from Figure 2(c). 

5 Conclusion 

To conclude, in this paper we have studied the equilibria and the phase structure of the 
asymptotically flat dilatonic black p-branes in a fixed cavity in arbitrary dimensions D 
in a grand canonical ensemble. For this, we have considered both the temperature and 
the potential fixed at the wall of the cavity and compared with our previous study of the 
same system in a canonical ensemble [1] for which the charge enclosed in the cavity, rather 
than the potential at the wall of the cavity, is fixed. We computed the Euclidean action 
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corresponding to the dilatonic black p-brane solution of .D-dimensional supergravity with 
maximal supersymmetry with proper consideration of the above fixed quantities at the 
wall of cavity in this ensemble. To the leading (or zeroth-loop) order this action is related 
to the grand potential or Gibbs free energy of the system and is an essential entity for 
the stability analysis. This action has two variables, namely, the horizon size r + (related 
to the variable x) and the charge Qd (related to the variable q). At the stationary point 
of this action, we derived the expected equations for which x and q have to satisfy with 
both the temperature and the potential fixed at the wall of cavity. The second derivatives 
of the action with respect to x and q at the stationary point can be used to analyze the 
local thermal stability of the system in this ensemble. 

Unlike in the canonical case, the directly obtained action here is a function of two 
variables and this makes the stability analysis a bit involved. However, there exists also 
an analog of one- variable analysis as in the canonical case and we have shown in the text 
that the condition so obtained for the stability of the equilibrium states as given in (39) 
is the same as that of the rather complicated two- variable case. For the two- variable 
analysis given in the Appendix, we have made use of a trick to use the 6, ip variables as 
defined in (30) instead of original x, q variables. This makes the analysis a bit simpler 
and straightforward. The stability condition can also be expressed as those given in (41). 
So, for a given potential satisfying the constraint given in (41), only the horizon size lying 
in the range also given in (41) gave locally stable system. 

We then obtained the phase structure of the black p-branes in this ensemble. We 
found that locally stable black brane phase exists only when the minimum temperature 
of the system given in (47) satisfies the constraint (48). In this case, below the minimum 
temperature, there is no stable black brane phase and only stable phase is the 'hot flat 
space'. Above this temperature and below the value T(l) given in (44), there are two 
black brane phases, the larger one is locally stable and the smaller one is unstable. The 
locally stable black brane phase becomes globally stable only above the temperature T g 
given in (55) and below T(l). Below T g and above T min (given in (47)), the stable black 
brane is only locally stable and will eventually make a transition to the 'hot flat space'. 
Finally, above T(l) there is no stable brane phase and the only stable phase is again the 
'hot flat space' phase. We also commented on the similarity of this phase structure with 
that of the zero charge canonical case except at one end of the x variable, namely, near 
x — 1. This structure is reminiscent of the Hawking-Page phase transition of the AdS, dS 
and flat black holes discussed in [8, 12, 13]. However, unlike in the canonical case, here we 
found that at a given temperature when the stable phase exists it does not come in pairs 
- only a single stable phase appears. So, there is no thermodynamical phase transition 
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between two stable black brane phases analogous to the van der Waals-Maxwell liquid-gas 
phase transition like that appeared in the canonical case. 
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Appendix: The two-variable stability analysis 



In this appendix we perform the more general yet more complicated two- variable analysis 
of the stability of black p-branes and show that the same local stability condition derived 
from the simple one-variable consideration as given in (39) continues to hold. 

It is clear from (30) that both b(x, q) and p(x, q) are smooth functions of x and q and 
so, in principle they can be inverted to give x and q as functions of b and tp. Although 
we will not need their explicit form, however, we will need the form of b as a function of 
x and tp. For that we will eliminate q from b(x, g)-equation given in (30) and express it 
as a function of x and tp instead. From the tp equation in (30) we obtain q 2 /x 2 as, 



.r- 



1 - (1 -ip 2 )x 



(64) 



Now substituting (64) in b equation in (30) we obtain 

b 



Xd [1 - (1 - p 2 )x} 2 



i i 

d 



The reduced action (31) can now be re-expressed as 



(65) 



Ie(x,<p) 



(d + 2)y + d(l- x(l - tptp)) y- 1 - 2(d+ 1) + bd(y 2 - l)y 



where we have used (64) and (65). In the above, we have also defined 

y = v 7 !- tp 2 ), 



(66) 



(67) 
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with x(b,ip) determined by (65). We will use both of (65) and (66) later. 

We now expand the reduced action (31) at the stationary point determined by the 
equation (29) with x = x and q = q to quadratic order as, 



lE(x,q) = lE(x,q) + I 



( Zi - Zi){Zj - Zj) + 



(6? 



Zi — Z j , Z j — Z j 



where Z\ = x, z-i — q with i, j — 1, 2 and we have used the stationary conditions for the 
first order terms. In the above 

d 2 I E (x,q) 



hj — 



dzjdzi 



(69) 



The local stability of the system is determined by whether the quadratic terms in the 
expansion are always positive definite. This can be easily understood if we diagonalize 
the matrix 1^ and demand that each of the two eigenvalues Ai and A2 is positive definite. 
Now since, 



AiA 2 = det Iij, and Ai + A 2 = I xx + /, 



(70) 



Ai, A 2 > implies that both I xx + I qq > and det 1^ > 0. Now we will rewrite these two 
conditions such that our analysis becomes simpler as 17 , 
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det L~ 



> 0, 



> 



(71) 
(72) 



In order to understand the meaning of the condition (71) we have to evaluate I qq and 
then set it to greater than zero. From the form of Ie in (31) we obtain the first condition 
(71) as, 



d 2 h 



1 qq 



dq 1 



x=x,q=q 



bdip 



q 1- 



d(l 



^)(i 



x) 



> 0. 



(73) 



Once this is satisfied the thermodynamic stability is completely determined by the con- 
dition (72). A direct evaluation of I qq /detlij from the expression (in (x, q) variables) of 



17 It is not difficult to see that the conditions given in (71) and (72) automatically implies det Iij > 
and I X x + Iqq > . To see this note that we first have to have I qq > (if I qq < 0, the above two conditions 
can never be satisfied), then (72) implies det > and so, det /j. 

=>■ Ixx - Iqx/Iqq > =^ L > Iqx/^qq > aild this in turn implies I x 



-I 2 

qx 



Iqq(Ixx Iqx/Iqq) > 



Iqq > 0. 
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Ie in (31) is complicated and tedious and we will make use of a trick in evaluating it. 
This will make the analysis of stability much more elegant and simpler. For this, we will 
instead evaluate the l.h.s. of (72) by first going to the (b, tp) variable and evaluate I a b, 
where 

Iab = ^PeT (74) 

with a, b = 1, 2 and £i = y?, £2 = ^ and then relate this matrix to the original matrix 1^ 
in (x, q) variables by chain rules as follows, 



lab 



hip hb 

dq dx 

dip dip 

dq dx 

db db 





' dq 


dx \ 


) ( 


d-p 


dv 


dq 


dx J 




v db 


db J 



(75) 



where 'T' denotes the transpose of a matrix. After that we invert this matrix relation to 
obtain, 



dq dx 

dip dip 

dq dx 

db db 



r 1 

1 ab 



dq dx 

dip dip 

dq dx J ' 

db db 



(76) 



Again inverting (76) we will evaluate iy/det hj in terms of I a b and other known functions 
given in (76). Now to evaluate we first calculate 



diE 

dp 



dh 



( dx\ 



dx 

[b - b(x,q)] f(x,q) 
- bd(<p - tp{x,q)) 



dlA fdq_ 
dq I \dtp 

/ X 

( dx 

, (d + 2)p(x,q) - 

+■ : — 3 ( b ~ ?)) 



dq 
dp 



(77) 



and 



diE 
db 



d! E 

dx 



( dx\ 



[b - b(x,q)] f(x,q) 
bd(p - (p(x,q)) 



dq 

dx 
db 



dq 
db 



(d + 2)<p(x,q) 



(b - b(x,q)) 



dq 
db 



{7i 
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where we have used the form of Ie given in (31) and the function f(x, q) is defined as, 



/(x,g) = (l-x)-* 1-^ 



d + 2 



Thus we have from eqs.(77) and (78) 
d 2 I E -~dq d 2 h 



bd— — , 



dip 2 dip d(pdb 

d 2 I E = - b ~dq cPIe 
dbdip db ' db 2 



d + 2 I 



1 - SI 



x z 



>0, (79) 



f( f g) dx I (d + 2)<p dq 
nx ' q % + 1-2! dip' 

X 

-f(r n) 9x I ^+2)Cp dq 

I[x ' q) db + \-t db" 



SO) 



where the second derivatives are evaluated at b = b and tp — tp. So, from the above 
equations (80) we obtain I ab as, 



Of 



(81) 



where we have used 



M db 



-f(x,q) 



dx (d + 2)ip dq 



dip 1 — si dip 

X 

From the form of lab given in (81) we can calculate its inverse as, 



I' 1 



mq)f h + m% -m 



det I ab 



db 
-Mil 



db 

bd^ 

atp 



with det I ab having the from (as can be seen from (81)) 

~ f dq dx dq dx 



det 1^ = fbd . 

\ db dip dip db 

where we have made use of (82). Now substituting (83) in (76) we obtain, 



1 dq 1 dx 

bd dip M d<p 

1 dx (d+2)<p dx 1 dx 



where we also made use of (82). Inverting (85) we get 



lij — det I{j 



(d+2)p dx _ Idx _±dx 

fbd(l-^ 9 V f db bd dip 

_l_dx 1 dq 

bd dtp bd dtp 



12) 



13) 



34) 



55) 



(86) 
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From (85) we also obtain, det J^ 1 = (det Iij) 
get, 



det Iab/U'bd) 2 - Now from (86) we finally 



1 </'/ 



(d + 2)<p dx 



det La 



fbdh 



dtp 



1 dx 



17) 



We can compute dx/db and dx/dtp at the stationary point from (65) and obtain, 



dx 2dx [1 - 5(1 - <p 2 )} dx 2 ^ 






db g 


2- (d + 2)5(1 -<^ 2 ) 


9( P (1 - ^ 2 ) 


2- (d + 2)5(1 -<^ 2 ) 



(88) 



Substituting (88) in (87) we have from the condition (72) 

{d + 2)(p 2 



<-qq 



2dx [1 - 5(1 - <p 2 )] 



det Iij 



fb 



(d + 2)x(l-(p 2 ) 



d(l -(^ 2 )(1 -x) 



> 0. (89) 



Thus for the stability of the system the above condition has to be satisfied. Now since 
the second factor in the square bracket is positive definite by the first condition written in 
(73) it is clear that for (89) to hold the denominator of the first factor has to be negative 
definite. This gives, as promised, precisely the same condition as given in (39). 

Now it can be checked that once (39) is satisfied the first condition written explicitly 
in (73) is automatically satisfied. To see this we look at the numerator of the second term 
in the square bracket in the expression of I qq given in (73), 

1 



1-|+§5(1-^ 2 ) 
d d 



d 



d-2(l-x(l-y 2 )) 



[1 - 5 (1 - <p 2 )] + -= [(d + 2)5 (1 - ip 2 ) - 2 



(90) 



and this is positive definite since the first term is positive definite and the second term 
is positive definite by (39). Therefore, I qq is positive definite. This shows that (39) is 
equivalent to both the stability conditions given in (71) and (72). 
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